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After imposing the Gauss law constraint as an initial condition upon the Hilbert space of the 
Nambu model, in all its generic realizations, we recover QED in the corresponding non-linear gauge 
A^A^ — n 2 M 2 . Our result is non-perturbative in the parameter M for n 2 ^ and can be extended 
to the n 2 — case. This shows that in the Nambu model, spontaneous Lorentz symmetry breaking 
dynamically generates gauge invariance, provided the Gauss law is imposed as an initial condition. 
In this way electrodynamics is recovered, with the photon being realized as the Nambu-Goldstone 
modes of the spontaneously broken symmetry, which finally turns out to be non-observable 

PACS numbers: ll.15.-q, 12.20.-m, 11.30.Cp 



I. INTRODUCTION 



The possible violation of Lorentz invariance has recently received a lot of attention, both from the experimental 
and the theoretical perspective. In this latter case, the interest lies mainly in connection with possible effects arising 
from the drastic modifications of space-time at distances of the order of Planck length, suggested by most of current 
quantum gravity approaches. Experiments and astrophysical observations put severe bounds upon the parameters 
Oj| describing such violations, which are still being improved. 

In this work we revisit a different and interesting idea related to spontaneous Lorentz symmetry breaking, which 

_C is that photons (and gravitons) could emerge as the corresponding Nambu-Goldstone bosons from that process. We 

focus on the case of electrodynamics. Normally the null mass of both particles is explained by invoking some type 

I ' of gauge invariance. Since this almost sacred principle has undoubtedly been fundamental in the development of 

£> physics, it is very interesting to explore the possibility that it could have a dynamical origin[l[. This idea goes back 

&\ • to the works of Nambu [2[ and Bjorken [3] together with many other contributions |4|. Recently it has been revived 

*/"") \ in references [5HL0|. 

One of the most explored approaches starts form a theory with a vector field i? M endowed with the standard elec- 
trodynamics kinetic term plus a potential designed to break the Lorentz symmetry via a non-zero vacuum expectation 
value < B^ >, which defines a preferred direction n^ in space-time. This potential also breaks gauge invariance. The 
original model was proposed by Nambu |2| and recently has been generalized to incorporate the so called bumblebee 
models [111 ]. In general, the subsequent symmetry breaking obtained from the non-zero minimum of the potential, 
splits the original four degrees of freedom B^ into three vectorial Nambu-Goldstone bosons A^ , to be identified with 
the photon satisfying the constraint A^A^ = n 2 M 2 , plus a massive scalar field c, which is assumed to be excited at 
very high energies. 

It is noteworthy to recall here that the constraint A^A* 1 = ±M 2 was originally proposed by Dirac as a way to derive 
the electromagnetic current from the additional excitations of the photon field, which now has ceased to be gauge 
deg rees of freedom, avoiding in this way the problems arising from considering point-like charges [12| ■ As reported in 
[13j, this theory requires the existence of an ether emerging from quantum fluctuations, but does not necessarily imply 
the violation of Lorentz invariance, due to an averaging process of the quantum states producing such fluctuations. 

In the framework of the Nambu model, calculations at the tree level [2| and to the one loop level [?} have been 
carried out, showing that the possible Lorentz violating effects are not present in the physical observables and that 
the results are completely consistent with the standard gauge invariant electrodynamics. This is certainly a surprising 
result, which can be truly appreciated from the complexity of the calculation in Ref.[7[ and that certainly deserves 
additional understanding. The above mentioned calculations have been performed using an expansion of the non- 
linear interaction up to order 1/M and 1/M 2 , respectively. The basic question that remains is to which extent the 
non-perturbative Nambu model describes the standard effects of full QED, which have been experimentally measured 
with astonishing precision. 

In this work we recover and extend to the space-like and light-like cases the results regarding the equivalence 
of the time-like version (n^n^ > 0) of the Nambu model and QED, given in Ref. [lCj. Our proposal is to start 
with the full analysis of these theories at the classical level and then move on to the quantum regime. Even at 
the classical level the problem is a difficult one because of the presence of the non-linear constraint A^A^ = n 2 M 2 , 
although non-linear gauges have been successfully used in the quantization of Yang- Mills theories [14 |. Despite this, 



and selecting an approach which does not require to completely fix the gauge, we show that the Nambu model, with 
the Gauss law as an additional constraint, imposed as an initial condition, has exactly the same canonical structure 
of the electromagnetic field. This point becomes evident when both canonical structures are expressed in terms of the 
electric E and potential A fields. Since the quantum dynamics is determined by the canonical structure, this result 
warrants that both theories, in their non-perturbative form, are equivalent at the quantum level, thus yielding the 
same physical results. 

The paper is organized as follows. In Section II we review the Nambu model, showing that its space-like version arises 
from the spontaneous Lorentz symmetry breaking (SLSB) of a bumblebee model. The remaining sectors of the Nambu 
model are considered as a natural generalization of the space-like case, without specifying the mechanism leading to 
the SLSB. In Section III we show that taking the phase space of the Nambu model and imposing the Gauss law 
constraint on it as an initial condition, the resulting canonical structure coincides with that of the electromagnetic 
field in the nonlinear gauge given by the particular SLSB considered. In Section IV we show that the separate 
procedures for dealing with the Nambu model, presented in Section III, can be given a unified description using a 
Lagrange multiplier to impose the resolution of the non-linear constraint in each case. Finally, Section V contains the 
conclusions and final comments. Appendix A summarizes the canonical version of QED, which we use as a benchmark 
to state the equivalence with the different realizations of the Nambu model plus the Gauss law. Appendix B includes 
a specific example of the construction of the potentials in the non-relativistic version of the non-linear gauge, for the 
case of a constant magnetic field. Finally, in Appendix C we discuss the dynamics of the full SLSB bumblebee model 
in its space-like version and show how electrodynamics is recovered in this case. 

II. SLSB AND THE NAMBU MODEL 

As we will show in the following, the Nambu model can be motivated as the low energy limit of the bumblebee 
model (the metric is rj^ = diag(l, —1, —1, —1)) 



±B lu ,B'»'-?L(B lt B»)-±(V> IV ' 2 



C{B lx ) = — A B^- — {B^)--{B tl B»f-J^, d^ = 0, (1) 



with 

B^u = d^B,, — d v B^ L , (2) 

after the 5*0(1, 3) symmetry is broken. Following the standard steps we obtain the broken symmetry Lagrangian 

C = -~F^F^(l + -) +L(i + 2L)F ltt ,(Ai t ff'u>) 

^ V °o/ ^o V Co J 

-\n 2 (fyw) (0"w) + A AM" (d„u) (duu) 

-\ (-2m 2 n 2 ) u 2 - Xaouj 3 - ^ - JM" (l + ~), (3) 

where we have successively set B^ = A^ u (x)n"a(x) = ^A^(x)a{x), a{x) = <r +u(x), and F^ v — 9M„ — <9M M . This 
parameterization requires 

B,B" = <rV, A^ = aln 2 . (4) 

Here n v is a constant vector characterizing the symmetry breaking vacuum according to 

n a time - like : 5(9(1,3) -> SO{3), (5) 

n a space -like : SO(l, 3) ->- 50(1,2). (6) 

Let us emphasize that the light-like case (n 2 = 0) is forbidden because then B^B^ = 0, so that there is no spontaneous 
symmetry breaking since the potential turns out to be zero. The Goldstone modes are contained in the field A 11 , 
which is restricted by the condition (Q|, so that it contains only three degrees of freedom. The field uj(x) describes 
the massive fluctuations of <j{x) around the minimum <7q of the potential, where 

M u = VW. (7) 



In this way the parameter m 2 in ([T]) must be positive. We also require A > 0, in such a way that the potential is 
bounded at infinity. The minimum of the potential has to be real and satisfies 



tZ^iZ 



<?o = \J — >0, (8) 

which imposes n 2 < 0. Let us further consider the case where the effective theory described by ([3]) is in the energy 
range (E « M u ) where the field oj(x) is not excited. In this situation we have 



C = —rF^FV _ j^ AfiA » = __ < o. (9) 



It is interesting to remark that the choice n 2 < in j3|) also produces the right sign normalization for the kinetic term 
of the field uj. 

In this way we have shown that the spontaneous Lorentz symmetry breaking of the bumblebee model |TJ) leads to 
the Nambu model [2| only in its space-like version. We consider the remaining sectors of the Nambu model as an 
extension of this result, without deriving them from the spontaneous symmetry breaking of a specific model. Thus 
we replace the second equation in @ by the requirement 

A^A* = M 2 n 2 , M 2 > 0, (10) 

where n 2 can now be positive, negative or zero. In each case we will show that n^ indeed describes the vacuum around 
which the symmetry has been broken by the Nambu-Goldstone modes A^. 

The standard analysis of the model is made by explicitly solving the constraint A^A^ — M 2 n 2 and subsequently 
substituting the result into the Lagrangian (0), thus providing a theory analogous to the standard non- linear sigma 
model [15]. As we will show in Section IV, an alternative compact way of dealing with the model is to introduce the 
constraint via a Lagrange multiplier. 

In the following we examine the conditions under which this theory is equivalent to standard electrodynamics, in 
such a way that A^ can be interpreted as the photon field in the non-linear gauge A^A^ = M 2 n 2 . The choices of n 2 
provide alternative forms in which the constraint A^A^ 1 — M 2 n 2 is to be solved, in such a way that the remaining 
symmetry is manifest. For the time-like case we have 

A) = VM2 + A 2 , (11) 

with A = (A 1 , A 2 , A 3 ) being three independent degrees of freedom, and where the variable Aq is explicitly invariant 
under rotations. 

The situation in the space-like case is 



An = ^M 2 + (A) 2 - (A^f - (A 2 ) 2 , (12) 

where the independent variables are (A , A 1 , A 2 ) and the variable A3 is explicitly invariant under the corresponding 
50(1,2) subgroup. 

In the case n 2 = we introduce the following parameterization [2] 

( B A 2 +A 2 \ A ( B A 2 +A 2 \ 
Ao = c[l + — + — l I , A 3 = c 1 + I , (13 

which provides the expansion parameter c for the chosen dependent fields, which is analogous to M in the remaining 
cases. Here the independent variables are B(x), A\(x) and ^(x). 

We identify the vacuum (^4 M ) of each sector of the theory by setting equal to zero all space-time dependent fields 
in A p. In this way we obtain 

n 2 > 0, (A^) = (M, 0,0,0), 
n 2 < 0, (^) o = (0,0,0,Af), 
n 2 = 0, (4A, = (c,0,0,c). (14) 



or 



III. ELECTRODYNAMICS AS A CONSTRAINED NAMBU MODEL 

The first hint upon how the Nambu model is related to electrodynamics arises from the Lagrangian equations of 
motion. The variation of §§§ produces 

SC = (d fl F^ -J V )5A V . (15) 

and the equations of motion for each case are obtained looking at the variation of the independent variables 

A- 
time — like case : 5Aq = — -SAi, (16) 

Aq 

An A 

space — like case : 5A 3 = — — 8 An — -^-SA a , (17) 

which produces 

(d k F k0 - J°) 4 1 + (d^ - J*) = 0, i= 1, 2, 3, (18) 

(d k F k0 - J°) + (d,F» 3 - J 3 ) ^- = 0, 

^3 

(d^ a - J a ) - (d^ 3 - J 3 ) ^ = 0, a = 1, 2. (19) 

^3 

respectively. Thus we see that in order to recover electrodynamics from the Nambu model we must further impose 
the Gauss law diF l ° — J° = as an additional constraint. This cannot be done at the Lagrangian level because such 
a constraint involves the velocity Aj. 

Recalling that the quantum dynamics is determined by the canonical structure, our general strategy to prove the 
equivalence between the quantized Nambu model plus the Gauss law and QED will be to perform the Hamiltonian 
analysis of the former and to compare with the corresponding formulation of electrodynamics, summarized in Appendix 
A. 

A. The time-like case 

We start from the Lagragian density ©, where we have explicitly solved and substituted the non-linear constraint 
in the form A = %/A 2 4- M 2 , 

CwiAi) = i ({dtAi-di (VA 2 +Af 2 )) 2 - (V x aA - JoVA 2 + Af 2 - J t A\ (20) 

The canonically conjugated momenta IT are 

IT = -j- = doAi - diVAP + A* = E h (21) 

oAi 

leading to the canonical Hamiltonian density % = EiAi — C 

U = ^E? + ^B 2 + Eidiy/M 2 + A 2 + WA 2 + Af 2 + J Z A\ (22) 



with the basic non-zero Poisson brackets 



The Hamiltonian equations of motion turn out to be 



{A(x,«), Bi(y,*)} = M 8 (x-y)- ( 23 ) 



Ek + (d k E k - J ) f* + a2 + [V x B]. - ,T = 0, (24) 



where 

E i = A i -d i y/M 2 +A*, B i = e ijk d j A k . (25) 



In other words, we have recovered the Lagrangian equations of motion (|18p for the Nambu model in this representation. 
The next step is to implement the Gauss law constraint via an arbitrary function x, so that the modified Hamiltonian 
density is 

H = \v 2 + hi 2 - (diEi - J ) (y/M 2 + A 2 + x ) + JiA\ (26) 

Redefining the arbitrary function 

6 = X + VM 2 + A 2 , (27) 

we finally obtain 

H = i (E 2 + B 2 ) - 9 (diEi - J ) + J t A\ (28) 

with the corresponding first order action 

S= J (fyfEiAi-ny (29) 

In this way, the Hamiltonian density adopts the form (|A13[) of the Appendix A. The Poisson bracket relations also 
coincide with those of (|A14[) . We have thus shown that the imposition of the Gauss law into the Nambu Hamiltonian 
density (j2"2"j) . which is performed in (j2l)]) . does in fact lead to classical electrodynamics. Let us remark that, as will 
be shown in Section IV, the dynamics of the Nambu model, for any choice of n 2 , preserves the Gauss law in time, so 
that it is sufficient to impose it as an initial condition. 

B. The space-like case 

Here we start from ©, but making the substitution A 3 = y ' M 2 + A\ — A a A a 

1 



£ N (A°,A a ) = -(-E 2 T + E 2 3 -B 2 )-J A°-J a A a -J 3 JM 2 +A 2 -A a A a , a = 1,2, (30) 



where the independent degrees of freedom are now A ,A a . Since A 3 is just a shorthand for \JM 2 + A 2 , — A 2 -,, we 
have 

A 3 = ^A -^A a , a = 1,2, (31) 

A3 A3 

which leads to 

E 3 = 4^0-4^-3^0, (32) 

A3 A 3 

E a = A a -d a A . (33) 

The canonically conjugated momenta are 

dL N A (A ■ A a . \ A 

H a = %f = E a E 3 ^ = A a d a A Q f^A Q ^A b d 3 A Q ) *± t (35) 

dA a A3 \A 3 A3 J A3 

with the non-zero Poisson brackets now being 

fA (x,i),n (y,0} = £ 3 (x-y), {A a (x,t), n fc (v,i)} = <U<S 3 (x-y). (36) 



Solving for the velocities in (j34|) and (|35|) yields 



T 



A o = ( ^no + ^- (A Q n a + A a d a A + A 3 d 3 A ) ) , (37) 



A a = H a + d a A + -^A U . (38) 



We can now write the electric field in terms of the momenta as follows 



E a =U a + Aj^, E 3 = A 3 ^. (39) 

Aq Aq 



The Hamiltonian density 

1 



H = Tl Q A Q + U a A a - ^ (E| + Ej - B 2 ) + J°A + J a A a + J 3 yJ NP + A 2 - A a A a (40) 



can be written as 

H = S ( A <^ + ^1) + Jn a n a + ^A a u a 

l A I A 

-A (d a (j^Aa) + d 3 (j^M) + daOa 

+ ^B 2 (A , A a ) + J°A + J a A a + J 3 ^M 2 +A 2 -A a A a . (41) 

To make contact with the standard form of the Hamiltonian in QED we rewrite 



Itt tt , n o „ „ 1^2 i /n 



U a n a + ^A a U a = -El--[^jA a A a ^ (42) 

using (|39p , and realize that the Gauss law adopts the following form in terms of the independent variables 

diEi = d a E a + d 3 E 3 = d a H a + d a (A a ^-\ + d 3 ( A ^J-\ (43) 

Substituting (g2J| and (gSJ) in (gT|) we obtain 

U = \e\ + h? a + ^B 2 (A , A a ) - A (8^ - J ) + J a A a + J 3 A 3 . (44) 

The above Hamiltonian density corresponds to the Nambu model and the resulting equations of motion coincide 
with those obtained in (jT9")l . Eq. (|4"4")l has just the form of the standard Hamiltonian density for electrodynamics, 
except that A 3 , E 3 are notations for the corresponding functions in terms of the dynamical coordinates and momenta. 
Notice that in this case Aq is not a Lagrange multiplier, so that we are still missing the Gauss law to recover the 
electrodynamics. 

Next, we follow the established procedure by imposing the Gauss law to the Nambu Hamiltonian (14"4"|) via the 
arbitrary function \- The result is 

U = l -El+ l -El + \&{A ,A a ) 

~{A +x)(diE t -J ) + J a A a + J 3 A 3 , (45) 

The redefinition 6 = Aq + \ leaves 

H = X - (E 2 + B 2 ) - 9 {8^ - Jo) + J a A a + J 3 A 3 , (46) 



which is precisely of the form (|A13[) . As emphasized at the end of Subsection A, it is enough to impose the Gauss 
law as an initial condition. Now we use the transformations from the original space- like variables A 0l A a , n ,II a to 
the set Ai,Ej (i,j = 1,2,3). 



A a = A a , A 3 = J M 2 + Al - A a A a , (47) 



n „ n , w , , 



E a = U a +A a -^, E 3 = -^JM 2 + A 2 -A a A a , (48) 

Aq Aq v 

together with the brackets (|36| for Aq, A a , Ho, II a , to verify that we recover the brackets among Ai, Ej given in (|A14j> . 
which define classical electrodynamics. In other words, we have shown that the change of variables ([47jl , (|48| is a 
canonical transformation. We also have 

n A + n a A a = E l A l , (49) 

as expected. 

As a final comment in this section we observe that the possible singularities envisaged in Section IV of Ref. [Tyj , 
arising in the space-like case due to the possibility of Aq — and which are not present in the time-like case, do in fact 
appear in our construction via the factor 1/Aq which occurs in many of the intermediate steps in this Section, and also 
in the corresponding equations in Section IV. Nevertheless, such singularity is absent both in the final Hamiltonian 
and in the final Dirac brackets. 

C. The light-like case 

In this case the degrees of freedom are the fields B and A a , a — 1, 2, according to (fl~3|) . To simplify the expressions 
we redefine c + B — >• B, and thus we get 

A = -L (4B 2 + A a A a ) , (50) 

A 3 = -L (4B 2 - A a A a ) . (51) 

From here on Aq and A 3 are just notation for the functions given by the above equations, and their time derivatives 
are 

A = i U 3 B + ^A a A a \ , (52) 

A 3 = i (a B - \A a AA . (53) 

In terms of these degrees of freedom, the components of the electric field are 

E 3 = A 3 - d 3 A = i (a B - ^A a A a - Bd 3 Ao) , (54) 

E a = A a -d a A . (55) 

As in the preceding cases, the starting point is the Lagrangian density ©, which takes the form 

C N {A\A a )= l -(Vl + El-&) 

- Jo { B + §)- J "{ B -§)- jaA - (56) 

The canonical momenta of A a and B can be written 

TT - 9C -F Aa F TT - dC - A °F <W\ 

and satisfy the non-zero fundamental Poisson brackets 

{B(x,t),n B (y,i)} = ( 5 3 (x-y), {A a (x,i) , II h (y, t)} = <W5 3 (x - y) . (58) 



From (l57l) , the electric field in terms of the momenta is 



and the velocities can be expressed 



E a = n Q + ^n B , 

ZAq 



A„ 



e 3 = -j-n B . 



ZA 
B = -J- [Bn B + Bd 3 A + ]-A a (n a + d a A ] 

An \ Z 



The Hamiltonian density 



H = U a A a + U B B 



E 2 + Ei 



B 2 



Jq[B 



3. 

4B 



J 6 [B 



4B 



J a A a 



takes now the form 



H 



-n 



B 



-n 



Aa n B n a + 1b 2 



2 a 2A Q "* ' 2A ~"~" ' 2 



d a Tl a + d 3 



U B B 

An 



"•1^ 



J° 



A + J 3 A 3 + J a A a . 



Introducing the electric and magnetic field according to (|59|). we have 



B 



A a . 



e 2 = — (u B y + (u a y + -^u B u Aa 

Aq An 



and the Gauss law becomes 



d % F M -J° = d 3 E 3 + d a E a - J° 



<M^n B ) + a a (n. a + ^n fl 



j° 



(59) 

(60) 
(61) 

(62) 



(63) 



(64) 



(65) 



These relations imply that the Hamiltonian can be expressed, using the electric and magnetic fields and the Gauss 
law, in the compact form 



H = i (E 2 + B 2 ) - (diF i0 - J°) A° + J a A a + J 3 A 3 , 



(66) 



which corresponds again to the Nambu model. As in the previous cases, imposing the Gauss law as an initial constraint, 
we obtain the Hamiltonian of the electromagnetic field in the nonlinear gauge A^A^ = 



H = ]- (E 2 + B 2 ) - (d t F lQ - J°) 9 + J t A l 



(67) 



where 9 is an arbitrary function, corresponding to the primary constraint (diF lQ — J°) . Similarly to the space-like 
case, in the intermediate steps of this construction it is possible for singularities to appear because of the factor 1/B. 
However, as in the space-like case, such singularities are not present in the final Hamiltonian density and in the Dirac 
brackets. Once more, starting from the definitions (150H51[) together with (J54H55J) in terms of the canonical variables, 
and using the brackets ([551) it is a direct matter to recover the Poisson algebra of electrodynamics (|A14p in terms of 
the fields As and Ei. 



IV. UNIFIED DISCUSSION OF THE NAMBU MODEL 



The alternative possibilities to realize the spontaneous symmetry breaking can be jointly discussed by incorporating 
the non-linear constraint via a Lagrange multiplier A 



C 



1 



2^2\ 



F^F^ - J^A" + \{A^ - M 2 n 



(68) 



which yields the equations of motion 



dy.F^ -J v + 2XA V = 0, (69) 

A^A* - M 2 n 2 = 0. (70) 



From Eq. (|68|) we obtain the relation 



do (diF™ J") = -d k (^ (ft** J )) , (71) 

which implies that the imposition of the Gauss law as an initial condition guarantees that diF M — J° = remains 
valid for arbitrary times. It is also important to emphasize that the choice X(x, t = 0) = as an initial condition leads, 
according to Eq. (|68[) . to the requirement that all four Maxwell equations must be satisfied at t = 0, in particular the 
Gauss law. 

There are several possibilities to work out X(x, t). If we take a particular value for the free index in Eq. (|69p , v = o, 
to solve the Lagrange multiplier, it yields 

* = -^(d»F^-r), (72) 

(no sum over a) and thus we obtain 

A a {d^ v - r) - (diF ia - ,r) A" = o, (73) 

A^A* - A/ V = 0. (74) 

For cr = these equations correspond to the time- like case (fT5)l . and for a = 3 they produce the space-like equations 
of motion (|19p . This shows that the Lagrange multiplier formulation of the Nambu model is equivalent to the 
direct-substitution approach developed in Section III. If we introduce the definitions 

Ei = A, - d l A , B % = e ljk d J A k 1 (75) 

and impose an additional constraint, the Gauss law V • E — J° = 0, the above equations become Maxwell equations in 
the nonlinear gauge ((74)) . Let us emphasize that (|73)) contains the two cases corresponding to (|TT|) and (J12I) , together 
with the light-like choice n 2 = 0. 

The quantum dynamics of this model is determined by the canonical formalism. For the Nambu model such 
a formulation can be obtained in a very compact way, directly from the Lagrangian density (|68p. using the Dirac 
approach for singular Lagrangians. The coordinate fields are A^ and A, and the definitions of their canonical momenta 
leads to 

n l = ^ = F°\ ^ = ^^0, tt a = ^=0. (76) 

dA i ' 8A° 8X 

From here we can make explicit the time derivative of the spatial components of the vector field 

A i = -K i + d i A°, (77) 

and two primary constraints emerge 

4>l = 7T , 4>2 = 7T A . (78) 

Thus, the extended Hamiltonian density results 

H = -\*>y + ~FijF ij + JiA* + (Jo - d^ 1 ) A Q - X (A^A" - M 2 n 2 ) + un + vtt x . (79) 

The consistency conditions of the primary constraints lead to the secondary ones 

h = {ttx,H} = ApA* - M 2 n 2 = 0, (80) 

04 = W H} = 2XA + (d^ - Jo) = 0. (81) 
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In turn, the equations of motion of these last constraints fix the arbitrary functions u and v. Thus we have the set of 
second class constraints <f>i, 4>2, <f>3, and </>4, which determine a subspace where the Hamiltonian density becomes 



1 i 1 



H = --TTiTT* + jFyF** + (Jo - fllTT*) A° + JiA*. (82) 



To make explicit the dynamics it is also necessary to construct the Dirac brackets that define the symplectic structure 
of the model. The constraint 4>3 admits several realizations. One of these expresses the time component of A^ in 
terms of the spatial ones (the time-like case) 



A = VM 2 - A*Ai. (83) 
Here the spatial components of A^ are the degrees of freedom, and the resulting Dirac brackets are 

{AJ(x,t),7r l (y,t)} D = V lS 3 (x-y). (84) 

Then, the canonical momenta conjugated to the A 3 fields are directly the ttj. Another realization is obtained when 
we take a spatial component as the dependent field (the space- like case), such as 

A 3 = (M 2 + A°A + A a A a ) 1/2 , (85) 
where a = 1, 2. Now the Dirac brackets among the degrees of freedom [A , A 1 , A 2 ) and the momenta ttj become 

{ A Q (x, t) , Tra (y, t)} D = - ) (M 2 + Al (x, t) + A a (x, t) A a (x, t)) 1/2 ^ 3 (x - y) , (86) 



{A (x, t) , n a (y, t)} D = -^^S 3 (x - y) , (87) 



A {x,t) 

{A a (x,t),7r 3 (y,t)} D = 0, 

{A a (x, t) , n b (y, t)} D = ^S 3 (x - y) . (89) 

The momenta ttj are not the canonical conjugate momenta of the degrees of freedom [A a , A 1 , A 2 ) , but from the above 
Dirac brackets we obtain 

{A 3 (x,t), 7 r 3 (y,i)} I) = {(Af 2 +A (x,i)Ao(x,t) + A a (x,i)A Q (x,i)) 1/2 , 7 r3(y,t)} c = 5 3 (x-y), (90) 

and thus the A 1 and ttj fields satisfy the same Dirac brackets as in the time-like case 

{^(x,t), 7 r J (y,i)} I) -r,i ( 5 3 (x-y). (91) 

From the first equation in (|76j) we have iii = — Ei in both cases. Imposing the Gauss law 

(diEi - J°) = 0, (92) 

as an initial constraint on the phase space of the Nambu model, we obtain the extended Hamiltonian density for the 
electromagnetic field 

H = - (E 2 + B 2 ) - J • A - 6 (V • E - J°) , (93) 

where is an arbitrary function, with the usual brackets between the potential and the electric field 

{A, (x, t) , Ei (y, t)} - S t] S 3 (x - y) . (94) 

The constrained phase space of the Hamiltonian for the free Nambu model includes subspaces where (|82|) becomes 
negative, because of the term A°diir l . In general the existence of these subspaces leads to instabilities and spoils 
unitarity. When the Gauss law is imposed as an additional constraint to obtain electrodynamics in the non linear 
gauge, this term disappears from the constrained Hamiltonian, which becomes positive definite and yields a unitary 
theory. 
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V. FINAL COMMENTS 

The Nambu model can be motivated by the spontaneous Lorentz symmetry breaking in the space-like sector of the 
bumblebee model given by Eq. ([1]). It is defined by the Lagrangian density of electrodynamics plus the non- linear 
constraint A^A^ = n 2 M 2 and, in this sense, it is similar to the non-linear sigma model. The standard way of dealing 
with the Nambu model is by solving the constraints in their different choices according to the character of n^, which 
now is arbitrary, and then substituting the solution into the Lagrangian. In this way, it is clear that the Nambu 
model contains three massless degrees of freedom, corresponding to the Nambu- Goldstone modes of the SLSB, and 
manifestly breaks gauge invariance. As shown by Nambu, at the classical and tree level Q, this model plus the Gauss 
law as an additional constraint, is equivalent to standard electrodynamics in the non-linear gauge A^A^ = n 2 M 2 . 
Further calculations including one loop effects 7] corroborate this conjecture. The above mentioned calculations are 
carried only up to order 1/M and 1/M 2 , respectively. On the other hand, Ref. [LOJ contains a discussion of the 
non-perturbative equivalence in the time-like case. 

Comparing with quantum electrodynamics a la Dirac, as described in Appendix A, where the Gauss law is imposed 
as a first class constraint upon the Hilbert space of the problem, we are able to show that the quantized Nambu model 
in the time-like and space-like cases, plus the quantized Gauss law imposed on the corresponding Hilbert space does 
indeed reproduce QED in a non-perturbative way with respect to the parameter M. This result is also extended to 
the light-like case (n 2 = 0). It is important to emphasize that the dynamics of the Nambu model guarantees the 
conservation of the Gauss law, in such a way that it is sufficient to impose it as an initial condition. This is the main 
result of this work. 

The proof goes as follows. For each case (time-like, space-like or light-like) we start from the corresponding 
canonical degrees of freedom, arising after solving the constraint A^A^ = n 2 M 2 , and obtain their Dirac brackets. 
After rewriting the standard electrodynamic variables Ai , Ei in terms of the canonical degrees of freedom and with 
the use of the related Dirac brackets, we prove that they satisfy the symplectic algebra of electrodynamics, given by 
the brackets (|A14[) . Moreover, each of the original Hamiltonian densities (f2"2"j) . (|41l) and (|63[) adopt the form (|A13[) of 
electrodynamics, after the Gauss law is imposed via an arbitrary function upon them. Since the quantum dynamics 
is determined by the canonical structure, we can now construct the corresponding quantum formalism. The fields 
Ai,Ei will satisfy the commutation relations of QED, arising from the brackets (IA14J1 . and the respective Hamiltonian 
densities will become the QED Hamiltonian after the Gauss law is imposed. Since this constraint is first class, we 
adopt here Dirac's point of view that the quantum theory is defined in the Hilbert subspace of physical states that 
are annihilated by the constraint. The fact that QED is recovered in the gauge corresponding to each solution of the 
non-linear constraint can be traced back to the particular initial form of the Nambu model used in each case. 

It is appropriate to remark that the imposition of the Gauss law in the otherwise non-gauge invariant Nambu model, 
is equivalent to restoring such invariance upon the full system, because this first class constraint is the generator of 
gauge transformations. For this reason a more accurate statement is that gauge invariance is indeed dynamically 
recovered from the Nambu model in order to reproduce QED, provided it is imposed as an initial condition. This 
statement coincides with the results in Refs. @, Q3- The results of Ref. [3] also support this point of view because 
the perturbative calculation performed there is carried on a Hilbert space where the photons are transverse to begin 
with, and therefore the Gauss law is imposed from the outset. 

As shown in Section IV, the separate procedures described in Section III can be given a unified description using a 
Lagrange multiplier to impose the resolution of the non-linear constraint in each case. In this way it is straightforward 
to show that for every realization of a spontaneous Lorentz symmetry breaking consistent with A^A^ = n 2 M 2 , and 
after imposing the Gauss law, we obtain electrodynamics in the corresponding non-linear gauge. Furthermore, using 
the Dirac method, we explicitly show that in the canonical formalism for the time-like and space-like cases, the 
canonical brackets among the electric and the potential field are the same as in electrodynamics. 

Nambu has also proposed that the solution of the associated Hamilton- Jacobi equation for a particle in a given 
electromagnetic field should provide, via the principal Hamilton Jacobi function, the function that allows to make the 
change of gauge from an initial one (here taken as the Coulomb gauge) to the non-linear gauge A^A^ — n 2 M 2 [2]. In 
Appendix B we have explicitly constructed an example of this procedure for the case of a constant magnetic field. We 
have considered the 1/M (with M — > oo) limit of the non-linear constraint in the time-like case, which is consistent 
with the non-relativistic limit of the Hamilton- Jacobi equation. 

Finally we discuss the full SLSB model, corresponding to the space-like case, arising from the bumblebee model 
((T|) , where the massive field becomes dynamical instead of being automatically frozen as is done in the Nambu model. 
The equations of motion again preserve the time evolution of the Gauss law. Its imposition as an initial condition 
produces a consistent solution for the massive field which leads to electrodynamics in the related space-like non-linear 
gauge. 
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Appendix A: Electrodynamics 

Let us briefly review the standard Hamiltonian formulation of electrodynamics. We start from the second order 
action 



S 2 (Ai,A ) = f d 3 yC 2 = fd 3 y ± ((A l -d l A^ 2 - (V x aA - J°A Q - J 1 A, 



(Al) 



The canonical momenta are 



dC 
IT = —r- = Ai-diAo = E h n = 0, 
dA t 



which produce the Hamiltonian density 



H =^n 2 + I (V x A) 2 - (diEi - J°) A Q + J 1 A, 



We also have the PB 



(A2) 



(A3) 



{A i (^,t),A j (y,t)} = 0, {Ei(x,t),Ei(y,t)} = 0, 

{i4 i (x,t),n i (y,t)} = %<5 3 (x- y ), {A (x,t),n ( y ,t)}-(5 3 (x- y ). 



(A4) 
(A5) 



The canonical theory is constructed via Dirac's method, due to the fact that the constraint E = IIo ~ is present. 
The extended Hamiltonian density is 



-He = ^E 2 + 1 (V x A) 2 -A (8.E, - J°) + J 1 A, + uU . 



(A6) 



The condition 



= E(x)= S(x), / d*yH E (v)>, 



(A7) 



leads to the Gauss law constraint 



Q = (diEi - J°) 



(A8) 



Finally we obtain Q — 0, in virtue of current conservation. In this way, the final first order action for electrodynamics 

is 



S 1 (A h E i ,A ) = d 3 yC 1 = d 3 y (EiAi-H), 



(A9) 



where 



H =-(E 2 + B 2 ) + J*Ai - A (8^ - J°) + uU a , 



(A10) 



with the two first class constraints 



rin ~ o, diEi - j° ~ o. 



(All) 



Normally one fixes 



n ~o, a ~o 



(A12) 



with an arbitrary function to be consistently determined after the remaining first class constraint is fixed, which 
yields 



n =^(e 2 + b 2 ) + fAi - e (diEi - j°) 



(A13) 
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The remaining PB are 

{A i (x,t),A j (y,t)}=0, {E i (x,t),E i (y,t)}=0, {A 4 (x,t),£,(y,t)} = %5 3 (x- y). (A14) 

The equations of motion are 

Ai = {Ai,H}, ^E = -A-V6. (A15) 

Taking Vx from the above equation we recover Faraday's law 

VxE = -VxA = B. (A16) 

The remaining equation arises from 

Bj = {Bi,JT},4E-VxB = J. (A17) 

Appendix B: Explicit gauge transformation in the non-relativistic case 

In general, it proves to be very difficult to go from a known gauge to the non-linear gauge A^A^ = M 2 . Let us 
assume that we start from the Coulomb gauge with potentials A^, then the required potentials in the non linear gauge 
are obtained by the transformation 

A^ = A^ + d^A, (Bl) 

together with the equation 

A^+d^y =M 2 , (B2) 

which determines the required gauge function A. 

As suggested by Nambu, A can be identified with the principal Hamilton- Jacobi function S corresponding to a 
particle with charge q and mass m moving in the corresponding electromagnetic fields. In fact, the related Hamilton- 
Jacobi equation is 

d^S-^y = m 2 c 2 , (B3) 

in such a way that the relation between the two problems is 

c c 2 m 
A = — S, M = . (B4) 

q q 

Here we give an example of such a construction for the time-like case, in the non-relativistic approximation. The 
non-relativistic Hamilton- Jacobi equation is 

(d t S-q<S>)=-^(d l S+-A l )\ (B5) 

where 

S = S T mcH. (B6) 

The minus sign choice in the square root arises from the corresponding form of the Hamilton- Jacobi equation for the 
free particle 

dtS+^(dSf=(). (B7) 

On the other hand, the M — )• oo limit of the non-linear gauge condition A^A^ — M 2 , is 

AotM = ±4l. (B8) 
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The constant M can be reabsorbed in Aq, so that Aq — > Aq =p M and the gauge condition is 



A^±§r f (B9) 



In this way, the non relativistic gauge function A is determined by 



- 1 \ U l -dAY „ 

* + -dt\) = - V 2M . ^ = (<D,A), (BIO) 



which is rewritten as 

2 



( g $ + £ ftA )=- ^ ° J ■ (Bit) 

Comparing (|B11[) with (IB5|) we obtain 

-2A-5, m=^, (B12) 

c c^ 

which reproduces the relativistic relation between A and S. 

Let us consider the simple case of a constant magnetic field B = £>k. In the Coulomb gauge we have 

I = 0, A = ^B x r. (B13) 

Following Nambu's suggestion we use the Hamilton- Jacobi approach to determine the gauge transformation function 
A. To this end we calculate the Hamilton- Jacobi principal function S. As is well known, this function corresponds 
to the action of the system evaluated at fixed end points. We calculate 

5(x ,i ;xi,ti)= ( dtL(x(t),*(t)), (B14) 

J to 

just by substituting the solution of the equations of motion into the Lagrangian and integrating. We take the end 
point conditions as 

a; (to) - 0, 2/(t = 0), z(t o )=0, 

x(h) = xi, y(tt) = 2/1, z{t\) = z\. (B15) 



The Lagrangian is 



L (x, x, t) = -x 2 + ^A ■ x = -x 2 + <L_ (atf - j/i , (B16) 

2 c 2 lc 



together with the equations of motion 



x-— y = 0, y + —x = 0, z = 0. (B17) 

mc mc 



The solutions are 



1 ( 0i\ , , x 1 ( 9i \ ■ 

x[t) = — — I x\ + yi cot — I (coswt — 1) + — I x\ cot — — y\ I sinwt, 

1 / 6»i \ 1 / 0! \ . 

J/(*) = o I ^ + Vl cot T J smu}t + o 1 Xl cot T ~ yi j ( cos w< ~~ 1 ^ ' 



(B18) 
(B19) 



z(t) = —t. (B20) 

ti 

with uj = qB /mc and 9\ — ujt\. Substituting them in the Lagrangian we obtain 

X (x(t), x(t)) = -^ + — (a; 2 + y 2 ) esc 2 (^J cos wt, (B21) 
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which produces the final result 

S(xo, to; xi.ii) = —^ + — (xj + yl) cot l—j . 

A direct calculation verifies that S satisfies the Hamilton- Jacobi equation: 

dS 



1 

2m 



dS_ qBy\ / dS qBx 
'dx + ~2cj + Vdy ~2cT 



as 
dz 



at 



o. 



From Eg. (|B12| ) the gauge function results 



llz 



B 



2 c t 
Hence, the potentials in the non-linear gauge are 



A ( x >*) = -o — T- - T (* + W ) cot 



/Be 



V2Af 



A n 



1 <9A 1 M 



D 2 



dt 2 c 2 t 2 8M 



+ j^( x2 + y 2 ) csc i ^-- l 



(Be 

\2M' 



a d\ B . 

A x = A x + — — = \ V + a; cot 1 

dx 2 V \2M 



/Be 



9A 



A y = A y + — = +—[ x - y cot 



B 

~2 



\2M 



dX _ M z 
oz c t 

It is straightforward to verify that indeed we have 



(B22) 



(B23) 



(B24) 



(B25) 



A = —A 2 . 
2M 



(B26) 



Appendix C: The full SLSB model 

In this Appendix we explore the relation between electrodynamics and the full SLSB model arising from the 
bumblebee model given in Eq. ([1]). We emphasize again that this Lagrangian describes a physical system only after 
the spontaneous symmetry breaking is realized. For example, the mass parameter in Eq. ([1]) has the wrong sign for 
describing a massive vector field B^ . As we have shown in Section II, a correct realization is obtained only in the 
case n 2 < 0. We take n 2 = — 1 and the corresponding Lagrangian reads 



C{A„u,&) = 



--F^F^ ( 1 



(TO 



1 
(JO 



1 



00 



F„„ (A"0"w) 



+ - (dpu) (O^lu) + ^A»A» (d^uj) (d u u) 



-~ (2m 2 ) u: 2 - Xaoco 3 - ~w 4 - J^A" (l + ^-) + 6 {A^A* + a 2 ) , 



(CI) 



where we have introduced the non-linear condition via the Lagrange multiplier 0. Here a 2 = m 2 /A, m 2 , A > 0. 

Next we examine the equations of motion in order to search for dynamical conditions that could reproduce electro- 
dynamics. In order to simplify the dynamics it is convenient to adopt the following parametrization of the Nambu- 
Goldstome modes A,. 



(^(l — [ 1 H I An — A^, (^LtU — (J^Uu 

co/ cr 



J [l^ J U 



Oi/C- f_i, 



(C2) 



where we also recall the notation a = <jq +uj. This is motivated precisely by the way in which the symmetry breaking 
was introduced in Section II. As expected we have the identity 



4 4 ctq °o 2 



2a 2 ' 



(C3) 
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which reproduces de first two lines on the right hand side of (IC1[) . The result is 

£(<?„, w.fi) = -\c^C^ - m 2 u 2 - X<j uj 3 - ^lu 4 - J^C» + ifi (c^ + (<r + w) 2 ) , (CM) 

with a redefinition of the Lagrange multiplier 0. The equations of motion are 

fyC"*" - J" + fiC = 0, (C5) 

C^C" + (cro + w) 2 - 0, (C6) 

-2m 2 uj - 3Xcr uj 2 - Xuj 3 + (cr + w) = 0, (C7) 
where the last one can be rewritten as 

(a +u)p- Xuj (w + 2a )] = 0. (C8) 
This equation is satisfied for arbitrary Q when u = — era- In this case we have 

C M = (1+^)^ = (C9) 

and the remaining equations consistently lead to J v = 0, so that the Goldstone field has no dynamics. The physically 
relevant configurations arise when u> 7^ — <7q. Then Eq. (|C8[) reads 



Q - Xuj (lj + 2a Q ) = 0. (C10) 

Following the discussion in subsection C-l in Ref. [lOj we could explore the possibility of setting the initial condition 
fl(x,t = 0) = 0. Notice that this would require to set the four Maxwell equations, including the Gauss law, as initial 
conditions, according to Eq. (|C5[) . We prefer to think the other way round and we impose (piC M — J°) = as an 
initial condition. Since Eq. (|C5|) has the same form as Eq. (l69l) . the dynamics guarantees that the Gauss law with 
respect to the field C lQ is valid for all times, which in turn leads to Q(x, t) = 0, from the zero component of Eq. (|C5|l . 
That is to say, we now have 

(diC i0 - J°)(x,t) = ->n(x,t)=0. (Cll) 

In this way we also recover Maxwell equations for C^ v (x,t) : d^C^ v — J v = 0. Nevertheless, let us observe that we 
are interested in the Gauss law [diF M — J°) = together with Maxwell equations for F M „. To this end we consider 
the equation of motion (|C10[) which yields the solutions 

w = 0, -► C li = A ll , d^F^ = r . (C12) 

u = -2a , ->■ C M = -A MI d^ v = -J v . (C13) 



In both situations the gauge fixing condition 



A^ +a 2 = (C14) 



is satisfied. The case u) = describes the electromagnetic field, in the non-linear gauge (|C14[) . coupled to the current 
J v , while uj = — 2a corresponds to an electromagnetic field in the same gauge, but coupled to the current — J". 

Summarizing, we have shown that the imposition of the Gauss law as an initial condition in the full spontaneously 
broken space-like bumblebee model leads to electrodynamics in the non-linear gauge A^A^ + 0$ = 0. The situation 
is analogous to the Nambu model, except that here we have obtained the condition uj(x.t) — dynamically as the 
consequence of imposing the Gauss as an initial condition, together with satisfying the corresponding equations of 
motion. 
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